Pointing to the minimum scatter: the generalized scaling relations for galaxy clusters

We introduce a generalized scaling law, M_tot = 10^K A^a B^b, to look for the minimum scatter in reconstructing the total mass of hydrodynamically simulated X-ray galaxy clusters, given gas mass M_gas, luminosity L and temperature T. We find a locus in the plane of the logarithmic slopes $a$ and $b$ of the scaling relations where the scatter in mass is minimized. This locus corresponds to b_M = -3/2 a_M +3/2 and b_L = -2 a_L +3/2 for A=M_gas and L, respectively, and B=T. Along these axes, all the known scaling relations can be identified (at different levels of scatter), plus a new one defined as M_tot ~ (LT)^(1/2). Simple formula to evaluate the expected evolution with redshift in the self-similar scenario are provided. In this scenario, no evolution of the scaling relations is predicted for the cases (b_M=0, a_M=1) and (b_L=7/2, a_L=-1), respectively. Once the single quantities are normalized to the average values of the sample under considerations, the normalizations K corresponding to the region with minimum scatter are very close to zero. The combination of these relations allows to reduce the number of free parameters of the fitting function that relates X-ray observables to the total mass and includes the self-similar redshift evolution.Comment: 6 pages, 3 figures. MNRAS in pres

Weighing simulated galaxy clusters using lensing and X-ray

We aim at investigating potential biases in lensing and X-ray methods to measure the cluster mass profiles. We do so by performing realistic simulations of lensing and X-ray observations that are subsequently analyzed using observational techniques. The resulting mass estimates are compared among them and with the input models. Three clusters obtained from state-of-the-art hydrodynamical simulations, each of which has been projected along three independent lines-of-sight, are used for this analysis. We find that strong lensing models can be trusted over a limited region around the cluster core. Extrapolating the strong lensing mass models to outside the Einstein ring can lead to significant biases in the mass estimates, if the BCG is not modeled properly for example. Weak lensing mass measurements can be largely affected by substructures, depending on the method implemented to convert the shear into a mass estimate. Using non-parametric methods which combine weak and strong lensing data, the projected masses within R200 can be constrained with a precision of ~10%. De-projection of lensing masses increases the scatter around the true masses by more than a factor of two due to cluster triaxiality. X-ray mass measurements have much smaller scatter (about a factor of two smaller than the lensing masses) but they are generally biased low by 5-20%. This bias is ascribable to bulk motions in the gas of our simulated clusters. Using the lensing and the X-ray masses as proxies for the true and the hydrostatic equilibrium masses of the simulated clusters and averaging over the cluster sample we are able to measure the lack of hydrostatic equilibrium in the systems we have investigated.Comment: 27 pages, 21 figures, accepted for publication on A&A. Version with full resolution images can be found at http://pico.bo.astro.it/~massimo/Public/Papers/massComp.pd

Comparing the temperatures of galaxy clusters from hydro-N-body simulations to Chandra and XMM-Newton observations

Theoretical studies of the physical processes guiding the formation and evolution of galaxies and galaxy clusters in the X-ray are mainly based on the results of numerical hydrodynamical N-body simulations, which in turn are often directly compared to X-ray observations. Although trivial in principle, these comparisons are not always simple. We demonstrate that the projected spectroscopic temperature of thermally complex clusters obtained from X-ray observations is always lower than the emission-weighed temperature, which is widely used in the analysis of numerical simulations. We show that this temperature bias is mainly related to the fact that the emission-weighted temperature does not reflect the actual spectral properties of the observed source. This has important implications for the study of thermal structures in clusters, especially when strong temperature gradients, like shock fronts, are present. Because of this bias, in real observations shock fronts appear much weaker than what is predicted by emission-weighted temperature maps, and may even not be detected. This may explain why, although numerical simulations predict that shock fronts are a quite common feature in clusters of galaxies, to date there are very few observations of objects in which they are clearly seen. To fix this problem we propose a new formula, the spectroscopic-like temperature function, and show that, for temperature larger than 3 keV, it approximates the spectroscopic temperature better than few per cent, making simulations more directly comparable to observations.Comment: Submitted for publication in MNRAS; 15 pages, 10 color figures and 13 BW figures,mn2e.cls. High resolution figures available here: http://people.roma2.infn.it/~mazzotta/preprints/mazzotta.pd

Cool Core Clusters from Cosmological Simulations

We present results obtained from a set of cosmological hydrodynamic simulations of galaxy clusters, aimed at comparing predictions with observational data on the diversity between cool-core (CC) and non-cool-core (NCC) clusters. Our simulations include the effects of stellar and AGN feedback and are based on an improved version of the smoothed particle hydrodynamics code GADGET-3, which ameliorates gas mixing and better captures gas-dynamical instabilities by including a suitable artificial thermal diffusion. In this Letter, we focus our analysis on the entropy profiles, the primary diagnostic we used to classify the degree of cool-coreness of clusters, and on the iron profiles. In keeping with observations, our simulated clusters display a variety of behaviors in entropy profiles: they range from steadily decreasing profiles at small radii, characteristic of cool-core systems, to nearly flat core isentropic profiles, characteristic of non-cool-core systems. Using observational criteria to distinguish between the two classes of objects, we find that they occur in similar proportions in both simulations and in observations. Furthermore, we also find that simulated cool-core clusters have profiles of iron abundance that are steeper than those of NCC clusters, which is also in agreement with observational results. We show that the capability of our simulations to generate a realistic cool-core structure in the cluster population is due to AGN feedback and artificial thermal diffusion: their combined action allows us to naturally distribute the energy extracted from super-massive black holes and to compensate for the radiative losses of low-entropy gas with short cooling time residing in the cluster core.Comment: 6 pages, 4 figures, accepted in ApJL, v2 contains some modifications on the text (results unchanged

Thermal Conduction in Simulated Galaxy Clusters

We study the formation of clusters of galaxies using high-resolution hydrodynamic cosmological simulations that include the effect of thermal conduction with an effective isotropic conductivity of 1/3 the classical Spitzer value. We find that, both for a hot ($T_{\rm ew}\simeq 12$ keV) and several cold ($T_{\rm ew}\simeq 2$ keV) galaxy clusters, the baryonic fraction converted into stars does not change significantly when thermal conduction is included. However, the temperature profiles are modified, particularly in our simulated hot system, where an extended isothermal core is readily formed. As a consequence of heat flowing from the inner regions of the cluster both to its outer parts and into its innermost resolved regions, the entropy profile is altered as well. This effect is almost negligible for the cold cluster, as expected based on the strong temperature dependence of the conductivity. Our results demonstrate that while thermal conduction can have a significant influence on the properties of the intra--cluster medium of rich galaxy clusters, it appears unlikely to provide by itself a solution for the overcooling problem in clusters, or to explain the current discrepancies between the observed and simulated properties of the intra--cluster medium.Comment: 4 Pages, 3 Figures, Submitted to ApJ-Letter

Spectroscopic-Like Temperature of Clusters of Galaxies and Cosmological Implications

The thermal properties of hydrodynamical simulations of galaxy clusters are usually compared to observations by relying on the emission-weighted temperature T_ew, instead of on the spectroscopic X-ray temperature T_spec, which is obtained by actual observational data. Here we show that, if the intra-cluster medium is thermally complex, T_ew fails at reproducing T_spec. We propose a new formula, the spectroscopic-like temperature, T_sl, which approximates T_spec better than a few per cent. By analyzing a set of hydrodynamical simulations of galaxy clusters, we also find that T_sl is lower than T_ew by 20-30 per cent. As a consequence, the normalization of the M-T relation from the simulations is larger than the observed one by about 50 per cent. If masses in simulated clusters are estimated by following the same assumptions of hydrostatic equilibrium and beta-model gas density profile, as often done for observed clusters, then the M-T relation decreases by about 40 per cent, and significantly reduces its scatter. Based on this result, we conclude that using the observed M-T relation to infer the amplitude of the power spectrum from the X--ray temperature function could bias low sigma_8 by 10-20 per cent. This may alleviate the tension between the value of sigma_8 inferred from the cluster number density and those from cosmic microwave background and large scale structure.Comment: 6 pages, 3 figures, to appear in the proceedings of the Rencontres du Vietnam "New Views on the Universe

Massive Halos in Millennium Gas Simulations: Multivariate Scaling Relations

The joint likelihood of observable cluster signals reflects the astrophysical evolution of the coupled baryonic and dark matter components in massive halos, and its knowledge will enhance cosmological parameter constraints in the coming era of large, multi-wavelength cluster surveys. We present a computational study of intrinsic covariance in cluster properties using halo populations derived from Millennium Gas Simulations (MGS). The MGS are re-simulations of the original 500 Mpc/h Millennium Simulation performed with gas dynamics under two different physical treatments: shock heating driven by gravity only (GO) and a second treatment with cooling and preheating (PH). We examine relationships among structural properties and observable X-ray and Sunyaev-Zel'dovich (SZ) signals for samples of thousands of halos with M_200 > 5 \times 10^{13} Msun/h and z < 2. While the X-ray scaling behavior of PH model halos at low-redshift offers a good match to local clusters, the model exhibits non-standard features testable with larger surveys, including weakly running slopes in hot gas observable--mass relations and ~10% departures from self-similar redshift evolution for 10^14 Msun/h halos at redshift z ~ 1. We find that the form of the joint likelihood of signal pairs is generally well-described by a multivariate, log-normal distribution, especially in the PH case which exhibits less halo substructure than the GO model. At fixed mass and epoch, joint deviations of signal pairs display mainly positive correlations, especially the thermal SZ effect paired with either hot gas fraction (r=0.88/0.69 for PH/GO at z=0) or X-ray temperature (r=0.62/0.83). We discuss halo mass selection by signal pairs, and find a minimum mass scatter of 4% in the \PH model by combining thermal SZ and gas fraction measurements.Comment: 19 pages, 14 figures, submitted to Ap

Cosmological hydrodynamical simulations of galaxy clusters: X-ray scaling relations and their evolution

We analyse cosmological hydrodynamical simulations of galaxy clusters to study the X-ray scaling relations between total masses and observable quantities such as X-ray luminosity, gas mass, X-ray temperature, and $Y_{X}$. Three sets of simulations are performed with an improved version of the smoothed particle hydrodynamics GADGET-3 code. These consider the following: non-radiative gas, star formation and stellar feedback, and the addition of feedback by active galactic nuclei (AGN). We select clusters with $M_{500} > 10^{14} M_{\odot} E(z)^{-1}$, mimicking the typical selection of Sunyaev-Zeldovich samples. This permits to have a mass range large enough to enable robust fitting of the relations even at $z \sim 2$. The results of the analysis show a general agreement with observations. The values of the slope of the mass-gas mass and mass-temperature relations at $z=2$ are 10 per cent lower with respect to $z=0$ due to the applied mass selection, in the former case, and to the effect of early merger in the latter. We investigate the impact of the slope variation on the study of the evolution of the normalization. We conclude that cosmological studies through scaling relations should be limited to the redshift range $z=0-1$, where we find that the slope, the scatter, and the covariance matrix of the relations are stable. The scaling between mass and $Y_X$ is confirmed to be the most robust relation, being almost independent of the gas physics. At higher redshifts, the scaling relations are sensitive to the inclusion of AGNs which influences low-mass systems. The detailed study of these objects will be crucial to evaluate the AGN effect on the ICM.Comment: 24 pages, 11 figures, 5 tables, replaced to match accepted versio